Dirty Paper Coding for the MIMO Cognitive Radio 

Channel with Imperfect CSIT 

Chinmay S. Vaze and Mahesh K. Varanasi 



Abstract — A Dirty Paper Coding (DPC) based trans- 
mission scheme for the Gaussian multiple-input multiple- 
output (MIMO) cognitive radio channel (CRC) is studied 
when there is imperfect and perfect channel knowledge at 
the transmitters (CSIT) and the receivers, respectively. In 
particular, the problem of optimizing the sum-rate of the 
MIMO CRC over the transmit covariance matrices is dealt 
with. Such an optimization, under the DPC-based trans- 
mission strategy, needs to be performed jointly with an op- 
timization over the inflation factor. To this end, first the 
problem of determination of inflation factor over the MIMO 
channel Y = H X X + H 2 S + Z with imperfect CSIT is investi- 
gated. For this problem, two iterative algorithms, which 
generalize the corresponding algorithms proposed for the 
channel Y = H[X + S) + Z, are developed. Later, the neces- 
sary conditions for maximizing the sum-rate of the MIMO 
CRC over the transmit covariances for a given choice of in- 
flation factor are derived. Using these necessary conditions 
and the algorithms for the determination of the inflation 
factor, an iterative, numerical algorithm for the joint op- 
timization is proposed. Some interesting observations are 
made from the numerical results obtained from the algo- 
rithm. Furthermore, the high-SNR sum-rate scaling factor 
achievable over the CRC with imperfect CSIT is obtained. 

Index Terms — Cognitive radio, dirty paper coding, infla- 
tion factor, covariance optimization. 

I. Introduction 

THE cognitive radio channel (CRC) was introduced in 
[I]. A cognitive radio is a device that can sense its 
, environment in real time and can accordingly adapt its 
■ transmission strategy. These are of current interest be- 
cause of the dramatically high spectral efficiency they can 
achieve pp. In pp, the authors introduced a more general 
cognitive protocol under which the CRC is an interference 
channel with degraded message sets [2]. 

The Gaussian multiple-input multiple-output (MIMO) 
interference channel consists of two transmitter-receiver 
pairs with each transmitter having a message for its paired 
receiver and the received signals are defined via equations 
Y\ = HuXi + H21X2 + Zi, Y2 — HizX± + H22X2 + Z%. 
Here, .- =1 are the fading channel matrices of dimen- 

sions rj x tjj the transmitted signals X\ ~ CA/"(0,Ex) and 
X2 ~ CA/"(0, £2) are subject to the power constraint^ of 
tr(Ei) < Pi and tr(S 2 ) < -P 2 ; and Z 1 ~ CAf(Q,I ri ) and 
Z2 ~ CJ\f(0,I r2 ) are the additive noises [2]. The Gaussian 
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1 Notation: For a square matrix A, tr(A), \A\, and rank(A) denote 
its trace, determinant, and rank, respectively. For any general matrix 
A, A* and A + denote the complex-conjugate transpose and pseudo- 
inverse of matrix A, respectively. vec(A) denotes the vector obtained 
by stacking the columns of A. I m is the m X m identity matrix. Eh 
denotes expectation over the random variable H. 



MIMO CRC is defined as the Gaussian MIMO interfer- 
ence channel in which the second transmitter (correspond- 
ing to signal X2) or the cognitive transmitter (CT) knows 
the message (or the codeword) of the first or the primary 
transmitter (corresponding to signal X\) non-causally pp. 

An achievable-rate region for the Gaussian MIMO CRC 
has been proposed in In their coding scheme, the CT, 
because of its non-causal knowledge, acts as a relay to aid 
the primary receiver and also transmits its own message to 
its paired (or cognitive) receiver. It employs dirty paper 
coding (DPC) |4J to cancel the interference at the cogni- 
tive receiver due to the signals intended for the primary 
receiver. With the assumption that the required channel 
matrices are known at the transmitters and the receivers, 
it is further shown in [3] that their achievable-rate region 
includes points corresponding to the sum-capacity of the 
MIMO CRC under certain conditions. Now, to achieve the 
sum-capacity, an optimization over the transmit covari- 
ances is required. This problem is studied in [5] where the 
authors propose the so-called adaptive sum-power iterative 
waterfilling algorithm which computes the sum-capacity 
and the optimal transmit covariances. 

Although the MIMO CRC is increasingly being stud- 
ied under the assumption of perfect transmitter-channel- 
knowledge (CSIT) [6], [3], [5], [7j, etc., not many papers 
[5] exist which deal with the practically important scenario 
of imperfect-CSIT CRC. We find it timely to consider the 
aforementioned problem of covariance optimization under 
imperfect CSIT. Towards this end, one needs to seek an- 
swers to the following two questions: 

1) DPC at the CT under imperfect CSIT: Since the 
channel seen by the cognitive receiver is of the form Y = 
H\X + H2S + Z (this will become more clear in Section 
IIII[) . where S is the interference known non-causally at the 
CT but not at its receiver, it is imperative to first study 
the problem of DPC over this channel when there is im- 
perfect CSIT of Hi and H2 . This problem is equivalent to 
the determination of the optimal inflation factor (see [4]) 
under imperfect CSIT. We studied a similar problem for 
the fading dirty paper channel (FDPC) Y = H(X + S) + Z 
in [9] and developed two iterative algorithms for determi- 
nation of inflation factor. These algorithms significantly 
improve the prior attempts mentioned therein. The same 
problem for the channel Y = H1X+H2S+Z, which we call 
the Generalized FDPC (G-FDPC), has been considered in 
[10] , but only in the special case of (all) single-antenna ter- 
minals, and a suboptimal solution is proposed. We study 
this important problem in Section |TTJ of this paper. 

2) Covariance optimization under imperfect CSIT: The 
problem of covariance optimization is considerably in- 
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volved, even under perfect CSIT. In [5], rather than using 
the achievable-rate region of [3J , the authors formulate the 
problem in terms of an outer-bound (obtained in [3J) to 
the capacity region which includes points corresponding to 
the sum-capacity of the MIMO CRC. This is a non-convex 
optimization problem and is converted into an equivalent 
convex-concave game using the 'MAC-BC transformations 
[llj . Then, for the resulting optimization, an iterative nu- 
merical algorithm is proposed. Unfortunately, the algo- 
rithm can not always guarantee the optimum solution. 

The use of the outer-bound or the 'MAC-BC transfor- 
mations is not possible under imperfect CSIT. Also, unlike 
the perfect-CSIT case (under which the interference can 
be assumed to be canceled perfectly by DPC), under im- 
perfect CSIT, an additional optimization over the inflation 
factor needs to be performed jointly with the transmit co- 
variances. Furthermore, the problem becomes more com- 
plicated because the sum-rate optimal solution need not 
necessarily have the power constraints satisfied with equal- 
ity (this point is detailed later). Thus the imperfect-CSIT 
version of this problem is also quite challenging. 

A slightly different (due to a constraint on the rate of 
primary user) version of the problem is considered in [8], 
[12] for the CRC with all single-antenna terminals. The 
authors of [5] , [T2] consider the amplify-and- forward strat- 
egy for relaying at the CT, and in this sense, their scheme 
is less general than the one studied here. 

II. DPC OVER THE G-FDPC 

Motivation to study this problem will become more clear 
in Section IlIII But, as noted before, this is an important 
step in the overall joint optimization. The G-FDPC is 
defined via equation Y = H\X + H 2 S + Z. Here, Hi and 
H 2 are the channel matrices of dimensions r xt x and r x t S) 
respectively; the transmitted signal X ~ CJV(0, Sj) has 
a power constraint of P; the interference S ~ CJ\f(0, Eg) 
is known non-causally at the transmitter but not at the 
receiver; Z ~ CM(0, Hz) is the additive noise; and A, S, 
and Z are independent. Assume perfect receiver channel 
knowledge but imperfect CSIT 0. Assume |£ X UE Z | > 0; 
let tr(E s ) = Q, tr(E z ) = N. Define SNR = ^ B Select 
the auxiliary random variable (see |13j for definition) as 
U = X + WS, i.e., Costa's choice g] extended to the MIMO 
case, where the t x x t s matrix W is the inflation factor 0. 
Similar to [9], we obtain the achievable rate as given by 



R 



maxEw log 
w 



|Ejc||E z + H{£, X H* + H 2 Y. S H 2 



with M = 
and H = 



E* + 



WT. S W* 
lH 2 T, s W* 



IMI 



V x Hl + WH S H* 



(1) 



Ez + HiY*xH\ + H 2 Y,sH 2 



Hi H 2 \- The above rate expression is valid only 
if I Ex I > 0. The case of Ex = can be handled as in 
[14] . We define the no-interference upper-bound i? n oS as 

2 We assume that the transmitter only knows the distribution H. 
The case of partial CSIT can be handled similarly. 

3 Note that N is total noise power. 

4 Matrix W is called the inflation factor so as to be consistent with 
the terminology introduced by Costa [4]. 




Fig. 1 

Achievable Rates vs. SNR: Rician Fading. 



the rate achievable over the G-FDPC in absence of inter- 
ference (i.e., when Q = 0) or R noS = E Hl log |Sz+ gJ| xg * 1 ■ 
The problem of determination of inflation factor, i.e., the 
maximization in ([1]) is equivalent to minn/E^ log IMI. 
As noted in [9] , this is a non-convex optimization problem, 
and it seems intractable to obtain a closed-form solution. 
It is possible however to generalize our algorithms in [9] 
developed for the FDPC to the G-FDPC. Due to lack of 
space, we discuss here the basic idea and omit the details. 

In the first algorithm, we minimize the objective function 
stepwise, i.e., at each step, we minimize over only one row 
of W, while treating all other rows as constants. Note 
that only the k th row and the k th column of matrix M 
depend on the k th row of W. Therefore, the minimization 
over one row of W (while treating other rows as constants) 
can be done analytically if the objective function is upper- 
bounded by moving the expectation inside the logarithm. 
Thus, one iteration of the algorithm consists of successive 
(stepwise) minimizations over all rows of W, and these 
iterations are repeated until a good choice is obtained. 

In the second algorithm, we solve for the stationary 
point of the objective function, i.e., solve an equation 
-JrrErH! h 2 ] log|M| = 0. Using the obtained necessary con- 
ditions, an iterative algorithm is proposed. 

Numerical Results: Here, 'lb' denotes the rate achiev- 
able using W = 0, i.e., by treating the interference as 
noise. R denotes the rate achievable using the algorithms. 
In Fig. [TJ we take Hi and H 2 to be independent with 
their elements ~ i.i.d. CJ\f(p,<r 2 ) with \/j,\ 2 + cr 2 = 1 and 
_ i., 1(1+3) 



(i=\H 



V2 



When /i = 0, the all-zero inflation factor per- 
forms almost as well as the inflation factor obtained using 
the algorithms. However, as \fi\ increases, the algorithms 
outperform the simple choice of W = 0. This type of obser- 
vation was also made in [10] in the case of SISO G-FDPC. 
It is generalized here to the MIMO case. In Fig. O the fad- 
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Fig. 2 

Achievable Rates vs. SNR: Hi and H2 Correlated. 



ing coefficients are correlated, i.e., H 1: H 2 ~ CAf(0, 1) with 
p = E(H*H 2 ). As p increases from to 0.7, the algorithms 
perform better than simply setting W = 0. 

A considerable difference between R and i? no s is seen in 
Fig. [TJ It should be noted that i? no s corresponds to the 
perfect interference-cancelation, while the curve R is for 
no CSIT. The gap between the two can be bridged with 
the availability of partial CSIT. Additionally, R no s is loose 
in the high-SNR regime because of the difference in the 
achievable scaling factors of R and i? n oS (see Theorem 1). 

Loosely speaking, it appears that for DPC to perform 
significantly better than the naive scheme of treating the 
interference as noise, it is necessary to have the matrix 
E(vec(Hi)*vec(H 2 )) 'non-zero', i.e., to have Hi and H 2 
'correlated'. The 'more' non-zero the above matrix is (or 
the 'more highly' Hi and H2 are correlated), the greater is 
the improvement. We believe this to be the fundamental 
nature of DPC over the G-FDPC under imperfect CSIT. 
Also see the discussion following Theorem 1. 

III. Optimization over the Transmit Covariances 

As per the coding scheme of [3], let X2 = X21 + A 22 
where the signal X21 corresponds to relaying and is corre- 
lated with Xi while X22 is the signal intended for the cogni- 



tive receiver. Let 



Xi 
X21 



CAM 0,E = 



Si 

v* 



V 
£21 



and 



X 22 ~ CA/"(0,S 22 ). Also let £ = 7\T* and £ 22 = T 2 T* 



for 

some Ti and T 2 ; and A 22 = T 2 X' 22 with X' 22 ~ CAf(0,I t2 ). 
The CT would choose the auxiliary random variable as 
Xi 
X 2 i 

and X 21 . Hence E 2 = S 22 + S 2 i. 

Now the channel between the CT-receiver pair is Y 2 = 

Xi ' 



U = X! 22 



W 



where X 22 



is independent of X\ 



H 22 X 2 



[H12 H 2 



Xn 



Z 2 which resembles the G- 



FDPC. Therefore, using the algorithms of Section [TTJ we 
can determine the inflation factor to be used at the CT 
once S (or T\) and S 22 (or T 2 ) are specified. This explains 
the reason to first study DPC over the G-FDPC. 

Denote Hi = [H n H 2 i], H 2 = [H 12 H 22 ], and H = 
[HI H 2 ]*. Then the achievable sum-rate R sum = R P + R C 
under no CSIT and perfect receiver channel knowledge is 
given by equation ([2|) at the bottom of the page. 

Since W depends on T\ and T 2 , we need to optimize 
Rsum jointly over Ti, T 2 , and W, as mentioned earlier. 
However, since W can be determined given T± and T 2 , 
let us first consider the optimization of R SU m over Ti and 
T 2 for a given value of W; later the algorithm for the 
joint optimization can be formulated. Let us consider: 
maxT u T 3 Rsum, subject to tr(Si) < P p and tr(S 2 i + S 22 ) < 
P c . This is a non-convex optimization problem. To obtain 
the necessary conditions, we form the lagrangian J; and 
set -jjjr — and J^r — 0. We omit the details of differenti- 
ation and directly state the necessary conditions, as given 
by equations ([3]) and ((4]) at the bottom of the page, where 
Ap 1 and A^T 1 are the lagrange multipliers. 

Algorithm for the Joint Optimization (Alg. 1): 



1. Start with some initial choices and T. 



For 



these choices, determine using the algorithms 

discussed in Section ITT1 
. At the n th iteration, 

> Determine the transmit covariances: to this end, we 



set T| 



in) 



andT 2 (n) =A c52 (Tr-^,T, 
quired expectations are evaluated numerically. Find 
lagrange multipliers so as to meet the power con- 



Aj, 1^ 







A c It 2 

(n-i) 



5H J 1 i L 2 
(n-1) ^(n-l)y 



,w (n - 1] ) 



The re- 



R sum = E s <^ log 



In + HiTiTfH* + H 2 {T 2 T' 2 H 21 \ 
|7 ri + H 2 iT2T 2 H 21 



dog 



tl 22 l 2 l2 tl 21 + "2-tl-i 1 n 2 



I t2 +WTiT*W* 
H 2 iT 2 + H 2 TiT^W* I r 



T 2 H 22 



WTiT*_H* 

IT rp rp* TT* 1 TT rp rp* 1 I • 



(2) 



A p 1 h 1 




A„ 1 I t2 



T 1= E R < HlN^Hl + H^N^Ii; 



[W* H* 2 ]D 2 



W 
Hi 



Ti=gi(Ti,T 2 ,W) 7 



\- x T2 = E S {H* 21 N^H 2 iT 2 - H* 21 D^H 2 iT 2 + H* 22 N 2 1 H 2 2T 2 - [0 H^]D^[J U T 2 *F* 2 ]*} 



g 2 {Ti,T 2 ,W), 



■ • where Ni = I ri + HiT x T*Hi + H 2 iT 2 T 2 *H^, D l = I ri + H 21 T 2 T^H^, 



N 2 =Ir 



H 22 T 2 T 2 H 22 + H 2 TiT* H 2 , and D 2 is the block-partitioned matrix in equation ((21). 



(3) 
(4) 



-1 



straints. 

• For Ti and obtained above, determine W^ n \ 
3. Repeat the above step until the increase in the achiev- 
able sum-rate is negligible. 
The statement above regarding the determination of 
lagrange multipliers warrants a discussion. Note, the 
power transmitted by either transmitter increases with 
A PjC . Therefore, the feasible region for the A's is of the form 
< X P ,c < A™^, where (\™ ax , \™ ax ) is a point at which 
both the power constraints are satisfied with equality. One 
can expect the optimal point to be (A" la:E , A™ ax ) at which 
both the transmitters operate with the maximum available 
power. However, since the signal intended for the primary 
receiver is an interference for the cognitive receiver and vice 
versa, the sum-rate need not necessarily be a nondecreas- 
ing function of either \ p or A c . Hence, the optimal point 
for A's, i.e., A° p c * can be any interior or boundary point of 
the above rectangular region. Note, the choice of A's dic- 
tates the covariance matrices, and therefore the inflation 
factor. Considering the fact that only an algorithmic solu- 
tion is available for the inflation factor and all the required 
expectations need to be evaluated numerically, the prob- 
lem of determination of optimal A's looks intractable. In 
the numerical examples, we consider a suboptimal solution 
of solving the power constraints as strict equalities. 

We have developed one more algorithm (Alg. 2) for the 
joint optimization which serves as a lower-bound on the 
rate achievable using Alg. 1. 

1. Assume that T\ = 0. Determine T 2 to maximize R c 
subject to tr(E 22 ) = -f (note the equality here). 

2. For given T 2 , determine T% to maximize R p under the 
constraints that tr(Ei) = P p and tr(E 2 i) = : 5 L - 

3. For given Ti, determine T2 and W to maximize R c 
under the constraint that tr(E 22 ) = -f- 

4. Repeat Steps 2 and 3 above until the increase in the 
achievable rate is negligible. 

Thus R p and R c are maximized here greedily over T± and 
(T 2 , W), respectively. For these maximizations, the algo- 
rithm of joint optimization developed in [TJ] is used. 

IV. High-SNR Analysis: Scaling Factor 

Theorem 1: G-FDPC : Assume that the ratio f is con- 
stant as P — > 00, and the fading processes are such that for 
any positive semi-definite A, rank([Bi H 2 ]A[Hi H2]*) = 
min(r, rank(j4)) with probability 1. The high-SNR scaling 
factor achievable over the no-CSIT G-FDPC using DPC 
is independent of the choice of W, as long as W is cho- 
sen such that the term log | Ex + VFEglV*! scales in the 
high-SNR regime as t x log SNR. 

Thus, the naive scheme of treating the interference 
as noise (i.e., W = 0) achieves the optimal scaling fac- 
tor, which is given by min(r, rank(Ex) + rank(Es)) — 
min(r,rank(£,g)). Also note that there is no loss of gen- 
erality in choosing W to satisfy the condition stated in 
Theorem 1 because Ws that do not satisfy this condition 
can achieve only a suboptimal scaling factor. 

The intuition detailed in the paragraph just preceding 
Section [TTI] can explain the result of Theorem 1. Con- 



sider the FDPC with Ex > 0, or equivalently, the G- 
FDPC with H 1 =H 2 = H and t x =t s = t. The high-SNR 
scaling factor of win(t,r), which is equal to that of the 
corresponding no-interference upper-bound, is achievable 
with the choice of W = It ['!• Next consider the FDPC 
with I Ex I = (let Ex = TT ); this channel is then equiv- 
alent to the G-FDPC Y = HiX' + H 2 S + Z with Hi = HT, 
H 2 = H, and X' ~ CM (0,1) (so Hi and B 2 not equal). In 
this case, the achievable scaling factor may not always be 
equal to that of the no-interference upper-bound; but in 
most cases, by making an appropriate choice for W (say, 
W — T + [H]), one can achieve a better scaling factor than 
that achievable with W = 0. Note, the G-FDPCs in the 
two cases above do not satisfy the assumption regarding 
Hi and H 2 made in Theorem 1. Finally, consider the G- 
FDPC that satisfies the assumption of Theorem 1 (for ex- 
ample, H\ and H 2 are independent and Rayleigh-faded) . 
Then, as per Theorem 1, there is no advantage in optimiz- 
ing over W as far as the scaling factor is concerned. Thus, 
the 'more highly' Hi and H 2 are correlated, the 'larger' is 
the increase in the scaling factor over that achievable by 
treating the interference as noise. 

Theorem 2: CRC : Assume that the channel matrices 
{Hij} are full rank and independent; and the ratio ^ 
remains constant. The high-SNR(= P p ) sum-rate scaling 
factor achievable over the no-CSIT CRC is given by 

j sum = max 7 P -I- 7 C with 

rank(S') , rank(S22) 

7 P = min(ri,rank(E')) — min(ri,rank(S 2 2)), 

7 C = min(r 2 ,rank(E')) — min(r 2 ,rank(E') — rank(E 22 )), 



where E' is the covariance matrix of 



Xi 

x 2 



and the max- 



imization is under the constraints of < rank(E') < ti +< 2 , 
and < rank(E 22 ) < min(< 2 ,rank(E')). 

Proof: We present the outline here and omit the de- 
tails. Given E and E 22 , R p achieves the scaling factor of 
7 P because the primary receiver treats the interference A 22 
as noise. 7 C is achieved by the choice of W = [0 T%] (see 
Theorem 1 above and Theorem 1 of [Hj). ■ 
The maximization in Theorem 2 is over only finitely many 
values; thus, can be done via exhaustive search. Note, the 
power constraints are not solved as strict equalities here. 

V. Numerical Results for the CRC 

In figures, 'ub' denotes the sum-rate achievable by op- 
timizing over the transmit covariances under the assump- 
tion that the interference is perfectly canceled at the cog- 
nitive receiver (i.e., the sum-rate with R c = Eh 22 log \I r . 2 + 
JT 22 E 22 /T2 2 |). We quantize each element of the fading ma- 
trices separately using an 'equally spaced level' quantizer 
as defined in [15]. In figures, if B — [By], then By denotes 
the number of feedback bits used per element of matrix 
H^. Further, By are independent. Alg. 1 is unfortunately 
sensitive to the initial choices. We take 4 to 5 initial choices 
in these examples and then select the best solution. 

In Fig. [3l we consider the CRC with elements of {B y } ~ 
i.i.d. A/"(0.6, 0.64). The improvement in the achievable 
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Fig. 3 

Achievable sum-rate vs. Pp. 




Fig. 4 

Achievable sum-rate vs. Pp. 



sum-rate with the introduction of partial CSIT is evident. 
Here, the scaling factor of 1 is achieved for R sum by let- 
ting the CT to use its entire power for relaying. Hence, the 
curve corresponding to B2 merges with that correspond- 
ing to the no CSIT at high SNR. For the CRC of Fig. 
IH we have the elements of {Hij} ~ i.i.d. Unif[0, 1]. For 
this CRC, as per Theorem 2, the optimal solution should 
achieve j sum = 2 with j p = 0. This fact can be easily seen 
from the plot. In Fig. [3 we have the CRC with elements 
of {Hij} ~ i.i.d. 7V(0, 1). It can be seen that Alg. 1 out- 
performs Alg. 2. However, in some cases, for example, the 
CRCs in Figs. [3] and 0] Alg. 2 does provide a relatively 
tight lower-bound. Coming back to Fig. [3] again, R sum 
achieves the scaling factor of 1 whereas according Theo- 
rem 2, j sum = 2. This is achieved by setting Si = 0, i.e., 
the primary transmitter needs to turn off its power. The 
apparent inconsistency here is because we have considered 
a suboptimal solution of solving the power constraints as 
strict equalities. This example emphasizes the importance 
of the problem of determination of A's. 




5 10 15 
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Fig. 5 

Comparison of Alg. 1 and Alg. 2. 



VI. Conclusion 

This paper is one of the earliest works that studies the 
imperfect-CSIT MIMO CRC. To the best of the authors' 
knowledge, it proposes for the first time a transmission 
strategy for the multi-antenna CRC with imperfect CSIT. 
En-route, brings into focus the problem of determination 
of A's. Furthermore, the paper derives an achievable high- 
SNR sum-rate scaling factor. It would be worthwhile to 
obtain the highest-achievable sum-rate scaling factor. This 
problem can be interesting; recall its counterpart for the 
Gaussian MIMO broadcast channel, a problem that is open 
even after serious attempts. More efforts are needed to 
answer these two open questions. 
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